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ABSTRACT 

In order to probe the order parameter symmetry of the heavy- fermion superconductor (HFS) CeCoIn 5 , we employ 
point-contact spectroscopy, where dynamic conductance spectra are taken from a nano-scale junction between a 
normal-metal (N) Au tip and a single crystal of CeCoIn 5 . The point-contact junction (PC J) is formed on a single 
crystal surface with two crystallographic orientations, (001) and (110). Our conductance spectra, reproducibly 
obtained over wide ranges of temperature, constitute the cleanest data sets ever reported for HFSs. The point 
contacts are shown to be in the Sharvin limit, ensuring spectroscopic nature of the measured data. A signature 
for the emerging heavy-fermion liquid is evidenced by the development of the asymmetry in the background 
conductance, starting at T* (~ 45 K) and increasing with decreasing temperature down to T c (2.3 K). Below T c , 
an enhancement of the sub-gap conductance arising from Andreev reflection is observed, with the magnitude of 
~ 13.3 % and ~ 11.8 % for the (001) and the (110) PCJ, respectively. These values are an order of magnitude 
smaller than those observed in conventional superconductors, but consistent with those in other HFSs. Our 
zero-bias conductance data for the (001) PCJ are best fit with the extended Blonder-Tinkham-Klapwijk model 
using the ci-wave order parameter. The fit to the full conductance curve of the (001) PCJ at 400 mK indicates 
the strong coupling nature (2A/kgT c = 4.64). However, our observed suppression of both the Andreev reflection 
signal and the energy gap indicates the failure of existing models. We provide possible directions for theoretical 
formulations of the electronic transport across an N/HFS interface in general, and the Au/CeCoIns interface in 
particular. Several qualitative features observed in the (110) PCJ provide the first clear spectroscopic evidence 
for the d x 2_ y 2 symmetry of the superconducting order parameter in CeCoIn5. 

Keywords: heavy-fermion superconductor, CeCoIns, Andreev reflection, Blonder-Tinkham-Klapwijk model, 
point-contact spectroscopy 

1. INTRODUCTION 

Probing the order parameter symmetry is of primary importance in investigating the pairing mechanism of 
unconventional superconductors. 1 Phase sensitive experiments such as Josephson junctions and SQUIDs 2 and 
single particle tunneling spectroscopy 3 have been playing important roles in manifesting the <i-wave symmetry 
of hole-doped cuprate high-T c superconductors (HTSs). Heavy-fermion superconductors (HFSs), another class 
of unconventional superconductors, have been studied mostly by point-contact spectroscopy 4-6 (PCS), partly 
because these materials are difficult to grow as thin films for tunnel junctions. 

Of the HFSs discovered, the relatively new family, CeMIn 5 (M = Co, Rh, Ir), have been of great research 
interest since they exhibit rich physical phenomena (Ref. 7 and references therein) such as quantum phase 
transitions 8, 9 and the Fulde-Ferrell-Larkin-Ovchinikov (FFLO) phase transition. 10, 11 There have been many 
reports on the symmetry of the superconducting order parameter in CeMIns. In particular, there exist substantial 
evidences for line nodes in the order parameter of CeCoIns. 12-15 The d- wave pairing symmetry is a most likely 
candidate but not conclusive yet since definite spectroscopic proof is still lacking. Goll et al. 16 and Rourke et al. 17 
have reported PCS data on CeCoIns, claiming that the superconducting order parameter has an unconventional 
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symmetry 16, 17 and even multiple components. 17 However, it can be shown that these data 17 may not reflect 
intrinsic properties of CeCoIn 5 . 18, 19 Furthermore, it remains controversial whether line nodes are located along 
the (lOO)-axis 13 or the (110) direction. 15 Addressing this issue is important for the clarification of either d x 2_ y 2- 
wave 13, 14 or d^-wave 15 symmetry. 

PCS, due to its simplicity and versatility as a spectroscopic tool, has been widely adopted for the investi- 
gation of both conventional and unconventional superconductors including HFSs. 4-6 In general, it can provide 
information on the density of states and the gap energy (A) through conductance data taken from a nano-scale 
junction between a normal-metal (N) and a superconductor (S). Blonder, Tinkham, and Klapwijk 20, 21 (BTK) 
formulated a theoretical model for the electronic transport across an N/S interface. The BTK theory provides 
clear descriptions on the transitional behavior from a metallic to a tunnel junction using the effective barrier 
strength, Z e ff, as a single parameter. Therefore, it has been playing a crucial role for the analysis of PCS data. 6 

If a quasi-particle (QP) is injected with energy lower than A from the normal-metal toward the superconductor 
in an N/S contact, it cannot enter the superconductor as a single particle since there are not available single 
particle states below A in the superconductor. This QP can transport to a superconductor by being retro-reflected 
as a quasi-hole and forming a Cooper pair with another electron. This Andreev reflection 22, 23 (AR) is a quantum 
mechanical scattering from a superconducting pair potential, conserving energy, momentum, spin, and charge. In 
a pure metallic N/S contact with Z e ff = 0, the zero- bias conductance (ZBC) is predicted to be doubled compared 
to that at high-bias (V 3> A/e), 20,21 so that the ZBC is enhanced by 100 %. A substantial enhancement of the 
ZBC due to AR was observed in point contacts containing conventional superconductors. 6,21,24 

According to the BTK theory, 21,24 the Fermi velocity mismatch acts as an effective barrier, thereby reducing 
the probability for AR. The effective barrier strength is given by Z e ff = [(1 — r) 2 / At+Zq] 1 / 2 , where r = vfn/vfs, 
the ratio of the Fermi velocities in the electrodes (note Z e f f remains invariant for r — > 1/r) , and Z a is the barrier 
strength due to an insulating layer. In mesoscopic semiconductor-superconductor junctions, results reported for 
Si-, GaAs-, InGaAs-, and InAs-based junctions with Nb counter-electrodes 25 could be accounted for using this 
formula. However, since it is not possible to separate the effects of an impurity- or disorder-induced barrier 
(Zo) at the interface from that of the Fermi surface mismatch in these systems, the accuracy of Z e ff remains 
inconclusive. An N/HFS point contact, due to the large disparity of the Fermi velocities, is expected to behave 
as a tunnel junction (Z e ff > 5). However, an enhancement of the sub-gap conductance (ESGC) due to AR has 
been commonly observed in many N/HFS point contacts, albeit suppressed in magnitude. 26-31 Deutscher and 
Nozieres 32 addressed this inconsistency between the theory and the experiments by proposing that the Fermi 
velocities entering in the ratio r are without a mass enhancement factor. 

Here, we report conductance spectra taken from point contacts on CeCoIns. We do observe an AR-induccd 
enhancement of the conductance but with heavily reduced magnitudes, implying that our results are consistent 
with reports on other HFSs. 26-31 Furthermore, the dependence of our data on the crystallographic orientations 
provides the first clear spectroscopic evidence for the <i-wave symmetry of the superconducting order parameter. 

2. EXPERIMENTS 

We have developed a PCS technique based on a combination of mechanical and piezoelectric mechanisms for 
making point contacts. 33 A schematic drawing of our PCS rig is displayed in Fig. 1. A finely polished Au tip is 
prepared by etching a Au wire in a concentrated hydrochloric acid with a dc pulse applied between it and a Pt 
counter-electrode. After the etching process, the surface is examined using an optical microscope and a scanning 
electron microscope (SEM). Only Au tips with smooth and clean surfaces are used. This is important to avoid 
the formation of non-ballistic and/or multiple contacts and degraded spectroscopic features. 

High quality CeCoIn 5 single crystals are grown using excess In flux. 34 Considering the tetragonal structure 
of CeCoIns, three kinds of point contacts are prepared with the surface having a different crystallographic 
orientation, named as (001), (110), and (100) point-contact junctions (PCJs). The (001) PCJs are formed on 
the largest surface of as-grown crystals since CeCoIn 5 single crystal is known to grow along the c-axis. Samples 
for (110) and (100) PCJs are prepared by embedding single crystals into a mold of a low temperature epoxy and 
cutting them such that the exposed surface is normal to the (110) and the (100) direction, respectively. These 
pieces are polished using alumina and/or diamond lapping films and silica colloidal suspensions down to 25 nm 
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Figure 1. A schematic drawing of the point-contact spectroscopy rig (not to scale). A point contact is formed by bringing 
a tip and a sample into contact using a fine screw. It is kept stable by adjusting the dc voltage to the piezoelectric bimorphs 
while being cooled down. The main body is made of Be-Cu to ensure a reversible movement when adjusted by the screw. 

particle size. The polished surface, examined using an optical microscope and an SEM, looks mirror-like shiny 
and smooth. The actual crystallographic orientation of the sample is checked by the X-ray diffraction. 

After a Au tip and a single crystal of CeCoIns are mounted onto the PCS rig, they are brought into contact 
by adjusting a fine screw at room temperature. Then, the PCS rig is attached to the 3 He cryostat. The point 
contact is kept stable by adjusting the dc driving voltage to the piezoelectric bimorphs during the cool-down 
process. The lowest achievable temperature is 300 mK and the maximum magnetic field is 12 Tesla. The dynamic 
conductance spectra of a PCJ are taken using the standard four probe lock-in technique as a function of bias 
voltage, temperature, and magnetic field. We report conductance spectra obtained from (001) PCJs and provide 
extensive analysis and discussion. 35 In addition, conductance data of (110) PCJs are presented with qualitative 
discussions. Results on (100) PCJs are not reported here, although some preliminary data have been obtained. 



3.1. Conductance of (001) Point-Contact Junction 

The dl/dV vs. V spectra obtained from a (001) PCJ are displayed in Fig. 2(a), normalized to the conductance 
at -2 mV. The PCJ was stable over a wide temperature range from 60 K to 400 mK. In the normal state, 
an asymmetry in the background conductance is seen to develop, starting at ~ 45 K. We attribute this to the 
emergence of a coherent heavy-fermion liquid. 36 This asymmetry is enhanced with decreasing temperature down 
to 2.6 K, below which it remains almost constant. This behavior is consistent with the observation that the 
relative weight of a coherent phase saturates below ~ 2 K. 36 These coincident behaviors of the background 
conductance with the emergence of the two fluids may provide important clues to our understanding of the 
electronic transport in CeCoIns. However, more detailed theoretical investigations are required. 

The constancy of the background conductance below T c enables us to perform quantitative analysis. That is, 
we can obtain symmetrized conductance spectra by factoring out the asymmetric part of the conductance data 
by that taken at 2.6 K, as shown in Fig. 2(b). Below T c , it is observed that the conductance near zero-bias is 
enhanced. This ESGC increases with decreasing temperature, showing a flat region at lower temperatures. We 
attribute this ESGC to AR at the Au/CeCoIns interface. At higher temperatures, this region becomes rounded 
due to the thermal broadening of the Fermi distribution function. As indicated in Fig. 2(b), the ZBC at the 
lowest temperature is enhanced by 13.3 %, which is much smaller than those in N/conventional superconductor 
point contacts with small Z e ff but consistent with the results for other N/HFS point contacts. 26,28 ' 31 We take 



3. RESULTS AND DISCUSSION 



this reduced ESGC as a ubiquitous penomenon in N/HFS point contacts, implying that there might be intrinsic 
mechanisms to cause AR to be suppressed at an N/HFS interface. 

In PCS measurements, it is important to ensure that conductance spectra are taken from a ballistic point 
contact to guarantee the spectroscopic nature of the data. Since it is not feasible to measure the size of a point 
contact directly, typically it is estimated using the measured resistance at a high bias (V » A/e), and Wexler's 
formula, 37 

+ <» 

where R is the resistance of the contact, K = l/a, a the contact radius, j(K) a smooth function of K, p and 
/ the resistivity and the electronic mean free path (EMFP) in electrodes. In the Sharvin 38 or ballistic limit 
(K > 1), ~/{K) -> 0.694, and R = R s = In the Maxwell 39 or diffusive limit (K -> 0), j(K) -> 1, and 

i?o = Rm = 2E- Since the low- temperature resistivity of Au is negligible compared to the reported value of 
CeCoIn 5 , 12 the contact size is estimated using only the latter. We also note that, for a point contact with non- 
zero Z e ff, Rq in Eq. (1) is related to the measured resistance, Rn, through the relation Rn = Rq(1 + Z e ff). 2A 
Rn ~ 1-1 and Z e ff ~ 0.365 (obtained through the analysis in a later section) are used. Since K is not known 
in advance, we estimate the contact size in both ballistic and diffusive limits, obtaining an upper limit for the 
contact size, 2a < 460 A. 

Movshovich et al. reported an estimation of the elastic EMFP to be 810 A at T c from the thermal conductivity 
data on CeCoIns. 12 Therefore, our measured point contact is much smaller than the EMFP at T = T c , indicating 
that our conductance spectra were taken from a point contact in the Sharvin or ballistic limit. An estimation 




Figure 2. Conductance vs. voltage data of a Au/CeCoInr, (001) point-contact junction. Curves are shifted vertically by 
0.05 for clarity, (a) Conductance data are normalized to the conductance at -2 mV. An asymmetry is seen to develop at ~ 
45 K, increasing with decreasing temperature, but remains constant below T c . (b) Conductance curves are symmetrized 
by dividing the data in (a) with that at 2.60 K. The enhancement of conductance near zero-bias, 13.3 % at 400 mK, is 
due to Andreev reflection. 



of the EMFP can be extended below T c using the same data 12 and the following thermodynamic relations for 
low-energy QPs in a d-wave superconductor. 40 



k/T cx p n r, p n cx T, t — l/vp, therefore, I cx k/T 2 , (2) 

where k is the thermal conductivity T the temperature, p n the density of normal QPs, r the QP lifetime, and 
vp the Fermi velocity. It is found that I increases exponentially with decreasing temperature, ranging 4-5 /xm 
at 400 mK, nearly two orders of magnitude larger than the contact size. We also estimate the inelastic EMFP 
based on the microwave conductivity data, 41 obtaining a lower limit of ~ 6500 A at 400 mK. Therefore, we can 
ensure that the measured contact is truly in the ballistic or Sharvin limit at low temperatures, even if we take 
into account the possibility of reduction of the EMFP in a point contact due to possible changes of the materials 
properties at the junction area from those of the bulk materials. These estimated values of the EMFP also show 
that CeCoIns is in the extreme clean limit (I ^> £, where £ is the coherence length), which together with the 
Pauli- limited upper critical field has been reported to be essential for the observation of the long-standing FFLO 
phase transition 10, 11 in this material. It is also clear that the arguments proposed by Gloos et al.f 2 attributing 
the suppressed ESGC to the non-ballistic nature of the contact, are not valid for our PCS measurements. 

As mentioned before, PCS data can provide information on the electronic density of states and the order 
parameter symmetry of a superconductor through analysis based on the BTK theory. Since the presence of 
line nodes in CeCoIns has been reported, 12-15 we choose the extended version of the BTK theory to a c?-wave 
superconductor by Tanaka and Kashiwaya, 43-45 hereafter called as the EBTK model. According to this model, 
the conductance is given by an integration of the conductance kernel over appropriate energy and momentum 
spaces. We have calculated the conductance of an N/d-wave superconductor junction as a function of bias 
voltage for three different cases, that is, the junction normal along the c-axis, lobe, and nodal directions of the 
superconductor. For a (001) PCJ, the conductance kernel is simplified due to its symmetry with respect to 
electronic trajectories. If the direction of QP momentum is strictly along the c-axis, which is the case of the 
extreme tunneling limit, no AR can occur. In an N/S point contact, typically Z e // is small. Therefore, if we follow 
the line of reasoning for the tunneling cone (the cone angle is inversely proportional to Z e jf), the conductance is 
obtained by the integration over the full half of the momentum space assuming that the momentum of the QP 
is distributed over this space: 

dl {v) = C Wfo* deTl dE d J^la s { El <t>) sin(2fl) 
iT^/c? d&rjv sin(20) 



where / is the Fermi distribution function, as(E,(f>) = ^jiqj-gg^j^p ^ > A = - ^\a\ CTjv = i+z i ' an< ^ 

Z = For the d-wave symmetry, A(T, (f) = A(T) cos 2(f>. In order to incorporate the effect of the QP lifetime 

broadening, we replace E = E' — iT, where r — Ti/t is the QP scattering rate, and take the real part of the 
kernel. 46 

Figure 3 shows conductance curves calculated for zero temperature and various values of Z e ff. Here, V is 
set to zero, a is the angle between the junction normal and the lobe direction of the <i-wave order parameter. 
In Fig. 3(a), it is seen that the conductance curves are identical for junctions with the normal along the c-axis 
and the lobe direction because of the reflection symmetry of the order parameter. The transitional behavior 
from AR to tunneling with increasing Z e ff is clearly demonstrated. AR do appear for a junction with small 
Z e ff, forming a conductance dip at the zero-bias with increasing Z e ff. On the other hand, a junction with the 
normal along the nodal direction exhibits strikingly different conductance features, as displayed in Fig. 3(b). 
For Z e ff = 0, the conductance curve looks the same as in Fig. 3(a), i.e., the usual AR-induced conductance 
is obtained. However, the conductance curve becomes narrower and sharper with increasing Z e ff, forming a 
zero-bias conductance peak (ZBCP) instead of a dip. The ZBCP in a tunnel junction made of a hole-doped HTS, 
a d-wave superconductor, has been frequently observed 3 and well understood theoretically. 47 It originates from 
the constructive interference between an incoming electron and an Andreev-reflected hole due to the phase shift 
of 7r in the order parameter at a surface of c?-wave superconductor whose normal is along the nodal direction. 
As a result, QPs form bound states at zero energy, called Andreev Bound States (ABSs). As demonstrated in 




Figure 3. Normalized conductance of a normal- metal/d- wave superconductor junction calculated using the extended 
BTK model with various effective barrier strengths, Z e ff. a is the angle between the junction normal and the lobe 
direction of the d-wave order parameter, (a) The junction normal is along the lobe or c-axis direction. With increasing 
Z e ff, a transitional behavior from Andreev reflection to tunneling is clearly seen, together with a dip in the zero-bias 
conductance for intermediate Z e ff. (b) The junction normal is along the nodal direction. The zero-bias conductance 
peak becomes narrower and stronger with increasing Z e f f . 

Fig. 3(b), the ABS-induced ZBCP is smeared out with decreasing Z e ff, 44 so that the conductance curve for 
Z e ff — is the same as in Fig. 3(a). 

Fitting to the data is performed by numerical integration of Eq. (3) using Z e ff, A, and T as parameters. 
Comparison of the experimental data at 400 mK in Fig. 2(b) to the calculation in Fig. 3(a) implies that Z e ff 
may have a value between 0.2 and 0.5. Although there exist many reports on line nodes in CeCoIns, 12-15 s-wave 
fitting is also performed for completeness. Since Rn remains almost constant below T c , it is reasonable to set 
Z e ff to a constant. The optimum value we obtain for Z e ff is 0.346. We also find that varying Z e ff as a 
function of temperature does not result in any better fits for either the s- or the (i-wave BTK model. The best 
fit curves using the s-wave BTK model are displayed as solid lines in Fig. 4(a). The flat region near zero-bias 
in the low temperature data are reproduced well by adjusting A and T. The calculated curve fits the data well 
near T c , whereas noticeable deviations are seen around the gap edge, these deviations increasing with decreasing 
temperature. This deviation is sometimes attributed to a local Joule heating effect due to the contact being in 
the thermal regime. 4 ' 48 However, this cannot be the origin in our point contact spectra since the contact is shown 
to be in the ballistic limit over the measured temperature range. The best fit values of A and V are plotted 
in Fig. 4(b) as a function of temperature. The gap energy extrapolated to zero temperature, A(0) = 404 fieV, 
gives rise to the ratio 2A(0)/fcsT c = 4.08, indicating the strong coupling nature in CeCoIn 5 in agreement with 
other experiments. 34 However, we note that T decreases with increasing temperature, which is unphysical and 
in contrast to usual observations 46,49 that T increases with increasing temperature. We attribute this behavior 
to the failure of the s-wave BTK model to account for our data. 

In the case of the (i-wave BTK model applied to a (001) PCJ, it is not possible to distinguish between d x 2_ y 2 
and d xy symmetry since the junction normal is along the c-axis. The best fit curve for the 400 mK data using 
the (i-wave model is displayed as a solid line in Fig. 5(a), together with the s-wave fit curve. Here, we point out 
that the shallow dip seen around -1.2 mV in the data is not an intrinsic feature indicative of the local heating 
effect, 4, 48 but an artifact caused in the normalization process due to an imperfect match of the background 
conductances at 400 mK and 2.6 K. The best fit values of the fitting parameters are, Z e ff — 0.365, T = 218 iieV, 
and A = 460 fieV, which gives the ratio 2A/fcsT c = 4.64, again implying the strong coupling nature. 34 As shown 
in Fig. 5(a), the (i-wave model gives a slightly better fit than the s-wave model, showing less deviation (albeit still 




Figure 4. (a) The best fit curves (solid lines) to the normalized conductance spectra (open squares) using the s-wave 
BTK model. Calculated curves and data are shifted for clarity. The temperature is, from the bottom to the top, 0.40, 
0.60, 0.80, 0.98, 1.12, 1.31, 1.52, 1.63, 1.75, 1.86, 1.97, 2.07, 2.15, 2.24, and 2.60 K, respectively, (b) Parameters used 
for the s-wave BTK fitting, A (filled squares) and F (open circles), together with the BCS energy gap (solid line) and 
r (dashed line) used in the fitting of the zero-bias conductance data in Fig. 5(b). A shows a reasonable temperature 
dependence, whereas F decreases with increasing temperature, which is unphysical. 

substantial) above the gap edge and reproducing a slight dip-peak feature near zero-bias. However, we cannot fit 
the conductance data over the whole temperature range without using T which shows an unphysical temperature 
dependence, i.e., decreasing with increasing temperature, similarly in the s-wave fit in Fig. 4(b). We interpret 
these results as a failure of the <i-wave EBTK model to explain the reduction in both the energy and the ESGC 
in our conductance curves. 

The ZBC would be least affected by any local heating effect. 4,48 Figure 5(b) shows that the ZBC vs. tem- 
perature data can be equally fit using both the s- and the d-wave models. The s-wave fit curve is obtained 
using A(0) = 349 [ieV and Z e ff = 0.346. However, again, the parameter T is required to have an unphys- 
ical temperature dependence, T(t) = 0.86A(0)(1 - i 3 /3), where t = T/T c , as plotted in Fig. 4(b). Com- 
bined together, these observations strongly indicate a breakdown of the s-wave BTK model in a (001) PCJ of 
CeCoIn 5 . This is not unexpected since most reports on the order parameter symmetry of this material sup- 
port the existence of line nodes. 12-15 Meanwhile, for the d-wave model, the fitting parameters are Z e ff = 0.365, 
A(T) = 2.35ksT c tanh(2.06y / T c /7 1 — 1), and T = 218 /J,eV, constant over the temperature range. This constancy 
of T in the d-w&ve fit is not unphysical, in contrast to the s-wave fit. Thus, we argue that the d-wave is a more 
likely pairing symmetry rather than the s-wave, consistent with the literature. 7, 12-15 

As seen above, however, the d-wave BTK model docs not fully account for our data taken over the full range of 
temperature below T c . We have investigated possible origins for this failure. First, we note that our observation 
of AR-likc conductance spectra is consistent with the PCS results on other N/HFS point contacts. 26-31 This is 
also in agreement with Deutscher and Nozieres' arguments, 32 in the sense that these N/HFS point contacts do 
not exhibit tunncling-like conductance features. In addition, the ESGC is commonly observed to be suppressed 
heavily by an order of magnitude, compared to that in N/conventional superconductor contacts. Our point 
contact spectra, although consistent with the data on other N/HFS contacts, constitute the cleanest data set 
over a wide temperature range. Since the ESGC is observed to accompany the superconducting transition, 
proximity 50 and pressure effects 48 ' 51 are ruled out as an origin of the suppressed AR. Likewise, the ballistic 
nature of the contact excludes the local heating effect 4,48 and the dominant Maxwell resistance. 42 As a result, 




Figure 5. (a) The best fit curve (solid line) to the conductance data (open circles) at 400 mK using the d-wave BTK 
model, together with the s-wave fit curve (dotted line). The d-wave fit gives rise to a slightly better result, reproducing the 
dip-peak structure near zero-bias, but the deviation is still substantial, (b) The zero-bias conductance vs. temperature 
data (filled squares) and the best fit curves using the s-wave (dotted line) and the d-wave (solid line) BTK models. Both 
fits are equally satisfactory, however, the temperature dependence of F is unphysical for the s-wave fit (see Fig. 4(b)), 
whereas acceptable for the d-wave fit (F = 218 ^eV). 

we claim that there must be intrinsic origins causing AR- induced conductance in N/HFS contacts to be reduced 
severely. 

Golubov and Tafuri 52 considered a breakdown of the Andreev approximation (retro-reflectivity) when A/Ep 
(Ep is the Fermi energy) is not negligible and/or the electrodes have layered structures as in the HTSs. Con- 
sidering the small Ep of 15 K for CeCoIns, 53 the non-retro-reflectivity may contribute to the reduced ESGC. 
Mortensen et a/. 54 considered mismatches in Fermi velocities and momenta and showed that the ESGC can be 
reduced in the case of large disparity. However, in their calculation, what is reduced is the SGC normalized by 
the normal state conductance instead of high-bias conductance (V ~> A/e). As shown in Fig. 5(a), the measured 
conductance above the gap edge is reduced compared to the calculated one. This may imply that the energy 
scale, which can be represented by A, is also reduced in CeCoIn 5 as well as the ESGC. Therefore, it is unlikely 
that the measured data can be accounted for with only mismatches in Fermi surface parameters. 52 ' 54 Looking at 
the two-fluid model, 36 one can formulate a model in which the effective barrier strength has different values for 
electrons and holes in each of the two fluids. 52 This kind of theory would naturally account for the asymmetry in 
the background conductance. On the other hand, it is found that conductance calculations, based on the EBTK 
models and taking into account spatial variations of the effective mass and the order parameter, just give rise 
to usual BTK conductances with proper scalings of the parameters. 55 Anders and Gloos 56 proposed that T in 
HFSs may be strongly dependent on the energy. Their calculated conductance curves for N/HFS point contacts 
based on the Green function formalism seems to be consistent with our data in the sense that both the energy 
gap and the SGC are reduced, compared to the BTK conductance for usual N/S point contacts. Nevertheless, 
this theory has a drawback that it is not easy to track down the physical mechanism clearly. More rigorous and 
detailed investigations of this model are necessary. One of other issues to be addressed is about the relevant 
time scale in regard to the QP broadening factor T: it has not been studied well for a multi-particle AR process, 
whereas plenty of investigations have been made for a single particle tunneling process. 46 The directionality of 
charge transport due to the quasi-two dimensional nature of the Fermi surfaces 57 and the effect of non-Fermi 
liquid nature of CeCoIn 5 8 should also be taken into consideration. 



3.2. Conductance of (110) Point-Contact Junction 

Figure 6(a) shows conductance spectra taken from a (110) PCJ. It is seen that the asymmetric background 
conductance remains almost the same over the temperature range from 2.58 K to 410 mK, which is also observed 
in the (001) PCJ. Therefore, the asymmetry in the background conductance seems to be an intrinsic feature of 
the Au/CeCoIns point contacts. Below T c , an enhancement of the conductance appears near zero-bias. Since 
Rn of this PCJ is <~ 4.7 CI, the contact size is about half of that for the (001) PCJ, ensuring the ballistic 
nature of the point contact. Similarly to the case of (001) PCJ, the conductance curves are symmetrized by 
dividing out the conductance data by the one at 2.58 K, as plotted in Fig. 6(b). The conductance curves at 
lower temperatures exhibit a change of the slope in the sub-gap region, which resembles the curve for small Z e f f 
in Fig. 3(b). It is remarkable that the cusp- like conductance feature near zero-bias persists up to the highest 
measured temperature below T c . 

The conductance data taken at the lowest temperatures are plotted together in Fig. 7. They are similar 
to each other in terms of the ZBC enhancement (13. 3 % vs. 11.8 %) and the energy scale for the conductance 
enhancement (<~ 1 mcV). However, the shape of the conductance curve is strikingly different near zero-bias: 
flat vs. cusp-like for the (001) and (110) PCJs, respectively. These different shapes are consistent with the 
calculated curves in Fig. 3, where clearly different behaviors of the conductance curve are observed with varying 
Z e ff. It is quite natural that the Au/CeCoIn 5 point contact may have non-zero Z e ff because of the disparate 
Fermi velocities, even after the Dcutscher-Nozieres theory 32 is taken into consideration. If the (110) PCJ were 
formed with a = 0, the conductance would then show a dip or flat structure near zero- bias, as in Fig. 3(a). The 
cusp-like structure in the conductance curve indicates that a ^ 0. The data are consistent with a = 7r/4, as 
shown in Fig. 3(b). This is consistent with the (110) PCJ being formed along the nodal direction of the order 
parameter, which supports the d x 2_ y 2-w&ve symmetry. 13 ' 14 As demonstrated in Fig. 3(b), the ZBCP of a PCJ 
with a = ir/4 increases rapidly with increasing Z e ff. The magnitude of the ZBCP is comparable to that of the 
(001) PCJ, implying that Z e ff is non-zero but small, maybe around 0.3, comparable to that for the (001) PCJ. 
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Figure 6. Conductance vs. voltage data of a Au/CeCoIns (110) point-contact junction. Curves are shifted vertically by 
0.05 for clarity, (a) Conductance data are normalized to the conductance at -5 mV. The asymmetry in the background 
conductance remains constant over the measured temperature range, similarly to the (001) point-contact junction, (b) 
Conductance curves are symmetrized by dividing out the data in (a) by the one at at 2.58 K. The conductance enhancement 
near zero-bias, 11.8 % at 410 mK, is due to Andreev Bound States. 
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Figure 7. Comparative plots of the symmetrized conductance data for the (001) and (110) point-contact junctions at 
400 mK and 410 mK, respectively. They are consistent in the sense that they have similar values for the zero-bias 
conductance and the energy scale for the conductance enhancement. The shape of the conductance curve in the sub-gap 
region is strikingly different from each other. The cusp-like shape at zero-bias and the slope change near ~ 400 fieV 
indicate that the (110) point-contact junction is normal to the nodal direction, consistent with the d x 2_ y 2 symmetry. 

The ZBCP in HFS junctions has been rarely observed, in contrast to HTS junctions. 3 This is partly because 
the techniques for thin film growth and junction fabrication are not well established. It is remarkable that our 
conductance spectra from the (110) PCJ show clear signatures of d x 2_ y 2-wave symmetry, thereby addressing the 
controversial issue 13, 15 on the location of line nodes. However, as in the HTSs, high quality tunnel junctions are 
desirable to provide more conclusive evidence for the ABS-induccd ZBCP. A complete and quantitative analysis 
on the conductance data from the (110) PCJ is deferred until a theoretical model is set up so that the mechanism 
for the suppressed AR, a common observation from both (001) and (110) PCJs, can be understood. 

4. CONCLUSIONS 

We have obtained conductance spectra from ballistic Au/CeCoIn 5 point-contact junctions with the normal 
along the c-axis and a6-plane directions, over wide temperature ranges. Our data constitute the cleanest sets 
ever reported for N/HFS point contacts. The data obtained from both junctions are consistent with each 
other in terms of the asymmetry in the background conductance, the magnitude of the ZBC, and the energy 
scale. Quantitative analyses of the conductance spectra for the (001) point-contact junction show that existing 
models cannot adequately describe the particle-hole Andreev conversion process at this interface. The gap 
energy extracted from the fit to the conductance curve at the lowest temperature implies the strong coupling 
nature. The temperature-dependence of a single point, the zero-bias conductance, is consistent with the d-wave 
order parameter symmetry. Both conclusions are consistent with the literature for CeCoIns. We propose that 
systematic corrections to the BTK model that go beyond the breakdown of the Andreev approximation and the 
re-normalized Fermi momenta may provide a framework for our future understanding of Andreev reflection at 
the N/HFS interface. Qualitative analysis of the conductance data for the (110) point-contact junction show 
clear signatures for the rf x 2_ y 2-wave pairing symmetry. Our data provide the first spectroscopic evidence for the 
order parameter symmetry and the orientation of line nodes in CeCoIn 5 . 

ACKNOWLEDGMENTS 



We are grateful to A. J. Leggett, D. Pines, V. Lukic, and J. Elenewski for fruitful discussions and to B. F. 
Wilken, A. N. Thaler, P. J. Hcntgcs, K. Parkinson, and W. L. Feldmann for experimental help. This work 
was supported by the U.S. Department of Energy Award No. DEFG02-91ER45439, through the Frederick Seitz 



Materials Research Laboratory and the Center for Microanalysis of Materials at the University of Illinois at 
Urbana-Champaign. 



REFERENCES 

1. K. H. Bennemann and J. B. Ketterson, eds., The Physics of Superconductors, Springer- Verlag, Germany, 
2003. 

2. D. J. Van Harlingen, "Phase-sensitive tests of the symmetry of the pairing state in the high-temperature 
superconductors — Evidence for d x i_ y i symmetry," Rev. Mod. Phys. 67, pp. 515-535, 1995; C. C. Tsuei and 
J. R. Kirtley, "Pairing symmetry in cuprate superconductors," ibid. 72, pp. 969-1016, 2000. 

3. M. Covington, M. Aprili, E. Paraoanu, L. H. Greene, F. Xu, J. Zhu, and C. A. Mirkin, "Observation of 
surface-induced broken time-reversal symmetry in YBa2Cu30y tunnel junctions," Phys. Rev. Lett. 79, pp. 
277-280, 1997. 

4. Yu. G. Naidyuk and I. K. Yanson, "Point-contact spectroscopy of heavy-fermion systems,"/. Phys.: Con- 
dens. Matter 10, pp. 8905-8938, 1998. 

5. H. v. Lohneysen, "Probing the energy gap of heavy-fermion superconductors," Physica B 218, pp. 148-156, 
1996. 

6. Yu. G. Naidyuk and I. K. Yanson, Point-Contact Spectroscopy, Springer Science+Business Media, Inc., New 
York, 2005. 

7. J. D. Thompson, M. Nicklas, A. Bianchi, R. Movshovich, A. Llobet, W. Bao, A. Malinowski, M. F. Hundley, 
N. O. Moreno, P. G. Pagliuso, J. L. Sarrao, S. Nakatsuji, Z. Fisk, R. Borth, E. Lengyel, N. Oeschler, G. 
Sparn, and F. Steglich, "Magnetism and unconventional superconductivity in Ce„M m In3„ + 2m heavy-fermion 
crystals," Physica B 329-333, pp. 446-449, 2003. 

8. V. A. Sidorov, M. Nicklas, P. G. Pagliuso, J. L. Sarrao, Y. Bang, A. V. Balatsky, and J. D. Thompson, 
"Superconductivity and quantum criticality in CeCoIns," Phys. Rev. Lett. 89, pp. 157004-1-4, 2002. 

9. J. Paglione, M. A. Tanatar, D. G. Hawthorn, E. Boaknin, R. W. Hill, F. Ronning, M. Sutherland, L. 
Taillefer, C. Petrovic, and P. C. Canficld, "Field-induced quantum critical point in CcColn 5 / Phys. Rev. 
Lett. 91, pp. 246405-1-4, 2003. 

10. R. A. Radovan, N. A. Fortune, T. P. Murphy, S. T. Hannahs, E. C. Palm, S. W. Tozcr, and D. Halll, 
"Magnetic enhancement of superconductivity from electron spin domains," Nature 425, pp. 51-55, 2003. 

11. A. Bianchi, R. Movshovich, C. Capan, P. G. Pagliuso, and J. L. Sarrao, "Possible Fulde-Ferrell-Larkin- 
Ovchinnikov superconducting state in CeCoIns," Phys. Rev. Lett. 91, pp. 187004-1-4, 2003. 

12. R. Movshovich, M. Jaime, J. D. Thompson, C. Petrovic, Z. Fisk, P. G. Pagliuso, and J. L. Sarrao, "Uncon- 
ventional superconductivity in Celrln 5 and CeCoIn 5 : specific heat and thermal conductivity studies," Phys. 
Rev. Lett. 86, pp. 5152-5155, 2001. 

13. K. Izawa, H. Yamaguchi, Y. Matsuda, H. Shishido, R. Settai, and Y. Onuki "Angular position of nodes 
in the superconducting gap of quasi-2D heavy-fermion superconductor CeColn 5 " Phys. Rev. Lett. 87, pp. 
057002-1-4, 2001. 

14. M. R. Eskildsen, C. D. Dcwhurst, B. W. Hoogenboom, C. Petrovic, and P. C. Canficld, "Hexagonal and 
square flux line lattices in CeCoIn 5 ," Phys. Rev. Lett. 90, pp. 187001-1-4, 2003. 

15. H. Aoki, T. Sakakibara, H. Shishido, R. Settai, Y. Onuki, P. Miranovic, and K. Machida, "Field-angle depen- 
dence of the zero-energy density of states in the unconventional heavy-fermion superconductor CeCoIns," J. 
Phys.: Condens. Matter 16, pp. L13-L19, 2004. 

16. G. Goll, H. v. Lohneysen, V. S. Zapf, E. D. Bauer, and M. B. Maple, "Hints at unconventional supercon- 
ductivity of the heavy-fermion superconductor CeCoIn 5 ,'Mcta Physica Polonica B 34, pp. 575-578, 2003. 

17. P. M. C. Rourke, M. A. Tanatar, C. S. Turcl, J. Bcrdeklis, C. Petrovic, and J. Y. T. Wei, "Spectroscopic 
evidence for multiple order parameter components in the heavy fermion superconductor CeCoIns," Phys. 
Rev. Lett. 94, pp. 107005-1-4, 2005. 

18. G. Sheet and P. Raychaudhuri, "Comment on spectroscopic evidence for multiple order parameters in the 
heavy fermion superconductor CeCoIn 5 ," http://arxiv.org/cond-mat/, pp. 0502632- 1-4, 2005. 

19. W. K. Park and L. H. Greene, "Comment on spectroscopic evidence for multiple order parameters in the 
heavy fermion superconductor CeCoIn 5 ," in preparation, 2005. 



20. G. E. Blonder, M. Tinkham, and T. M. Klapwijk, "Transitions from metallic to tunneling regimes in 
superconducting microconstrictions: Excess current, charge imbalance, and supercurrent conversion," Phys. 
Rev. B 25, pp. 4515-4532, 1982. 

21. G. E. Blonder, The I-V Curves of Superconducting Microconstrictions, Ph.D. thesis, Harvard University, 
1982. 

22. A. F. Andreev, "The thermal conductivity of the intermediate state in superconductors," 5W Phys. JETP 
19, pp. 1228-1231, 1964. 

23. G. Deutscher, "Andreev-Saint- James reflections: A probe of cuprate superconductors," Rev. Mod. Phys. 77, 
pp. 109 -135, 2005. 

24. G. E. Blonder and M. Tinkham, "Metallic to tunneling transition in Cu-Nb point contacts," Phys. Rev. B. 
27, pp. 112-118, 1983. 

25. W. M. van Huffclcn, T. M. Klapwijk, D. R. Hcslinga, M. J. de Boer, and N. van der Post, "Carrier transport 
in mesoscopic silicon-coupled superconducting junctions, "Phys. Rev. B47, pp. 5170-5189, 1993; F. Giazotto, 
P. Pinguc, and F. Bcltram, "Coherent transport in Nb/<5-Dopcd-GaAs hybrid microstructures," Mod. Phys. 
Lett. 17, pp. 955-971, 2003; A. Kastalsky, A. W. Kleinsasser, L. H. Greene, R. Bhat, F. P. Milliken, and 
J. P. Harbison, "Observation of pair currents in superconductor-semiconductor contacts," Phys. Rev. Lett. 
67, pp. 3026-3029, 1991; C. Nguyen, H. Kroemer, and E. L. Hu, "Contact resistance of superconductor- 
semiconductor interfaces: The case of Nb-InAs/AlSb quantum-well structures," Appl. Phys. Lett. 65, pp. 
103-105, 1994. 

26. G. Goll, H. v. Lohneysen, I. K. Yanson, and L. Taillcfcr, "Anisotropy in point-contact spectra of the heavy- 
fermion superconductor UP 1 3 ," Phys. Rev. Lett. 70, pp. 2008-2011, 1993. 

27. Y. DeWilde, J. Hcil, A. G. M. Janscn, P. Wydcr, R. Dcltour, W. Assmus, A. Menovsky, W. Sun, and L. 
Taillefer, "Andreev reflections on heavy-fermion superconductors," Phys. Rev. Lett. 72, pp. 2278-2281, 1994. 

28. G. Goll, C. Bruder, and H. v. Lohneysen, "Andreev scattering in the anisotropic heavy-fermion supercon- 
ductor UPt 3 ,"P%s. Rev. B 52, pp. 6801-6807, 1995. 

29. Y. DeWilde, T. M. Klapwijk, A. G. M. Jansen, J. Heil, and P. Wyder, "Quasi-particle lifetime broadening 
in normal-superconductor junctions with UPt3," Physica B 218, pp. 165-168, 1996. 

30. Yu. G. Naidyuk, H. v. Lohneysen, G. Goll, C. Paschkc, I. K. Yanson, and A. A. Menovsky, "Andreev 
reflection in point contacts between the heavy-fermion superconductor UPt3 and ordinary superconduc- 
tors," Physica B 218, pp. 161-164, 1996. 

31. C. Obermair, G. Goll, H. v. Lohneysen, I. K. Yanson, and L. Taillefer, "Surface-dependent Andreev reflection 
on UPt 3 ,"P%s. Rev. B 57, pp. 7506-7509, 1998. 

32. G. Deutscher and P. Nozieres, "Cancellation of quasiparticle mass enhancement in the conductance of point 
contacts," Phys. Rev. B 50, pp. 13557-13562, 1994. 

33. W. K. Park and L. H. Greene, "Construction of Cantilever-Andreev- Tunneling rig and application to su- 
perconductors," in preparation, 2005. 

34. C. Petrovic, P. G. Pagliuso, M. F. Hundley, R. Movshovich, J. L. Sarrao, J. D. Thompson, Z. Fisk, and 
P. Monthoux, "Heavy-fermion superconductivity in CeCoIn 5 at 2.3 K," J. Phys.: Condens. Matter 13, pp. 
L337-L342, 2001. 

35. W. K. Park, L. H. Greene, J. L. Sarrao, and J. D. Thompson, "Point-contact spectroscopy of the heavy- 
fermion superconductor CeCoIn 5 ," http://arxiv.org/cond-mat/, pp. 0409090-1-4, 2004. 

36. S. Nakatsuji, D. Pines, and Z. Fisk, "Two fluid description of the Kondo lattice," Phys. Rev. Lett. 92, pp. 
016401-1-4, 2004. 

37. G. Wcxlcr, "The size effect and the non-local Boltzmann transport equation in orifice and disk geome- 
try," Proc. Phys. Soc. 89, pp. 927-941, 1966. 

38. Yu. V. Sharvin, "A possible method for studying Fermi surfaces," Sov. Phys. JETP 21, pp. 655-656, 1965. 

39. R. Holm, Electric Contacts, Springer- Verlag Inc., New York, 1967. 

40. N. E. Hussey, "Low-energy quasiparticles in high-T c cuprates," ^4cfo. Phys. 51, pp. 1685-1771, 2002. 

41. R. J. Ormeno, A. Sibley, C. E. Gough, S. Sebastian, and I. R. Fisher, "Microwave conductivity and pene- 
tration depth in the heavy fermion superconductor CeCoIn 5 ," Phys. Rev. Lett. 88, pp. 047005-1-4, 2002. 



42. K. Gloos, F. B. Anders, B. Buschinger, C. Geibel, K. Heuser, F. Jahrling, J. S. Kim, R. Klemens, R. 
Miiller-Reisener, C. Schank, and G. R. Stewart, "Scaling behavior of point contacts between a tungsten tip 
and the heavy-fermion superconductors," J. Low Temp. Phys. 105, pp. 37-65, 1996. 

43. Y. Tanaka and S. Kashiwaya, "Theory of tunneling spectroscopy of cZ-wave superconductors," Phys. Rev. 
Lett. 74, pp. 3451-3454, 1995. 

44. S. Kashiwaya, Y. Tanaka, M. Koyanagi, and K. Kajimura, "Theory for tunneling spectroscopy of anisotropic 
superconductors," Phys. Rev. B 53, pp. 2667-2676, 1996. 

45. S. Kashiwaya and Y. Tanaka, "Tunneling effects on surface bound states in unconventional superconduc- 
tors," .Rep. Prog. Phys. 63, pp. 1641-1724, 2000. 

46. R. C. Dynes, V. Narayanamurti, and J. P. Garno, "Direct measurement of quasiparticle-lifetime in a strong- 
coupled superconductor," Phys. Rev. Lett. 41, pp. 1509-1512, 1978. 

47. M. Fogelstrom, D. Rainer, and J. A. Sauls, "Tunneling into current-carrying surface states of high-T c 
superconductors," Phys. Rev. Lett. 79, pp. 281-284, 1997. 

48. K. Gloos, "Andreev reflection at point contacts with heavy- Fermion UBei3?" P/m/s. Rev. Lett. 85, p. 5257, 
2000. 

49. Y. Bugoslavsky, Y. Miyoshi, G. K. Perkins, A. V. Bcrcnov, Z. Lockman, J. L. MacManus-Driscoll, L. F. 
Cohen, A. D. Caplin, H. Y. Zhai, M. P. Paranthaman, H. M. Christen, and M. Blamire, "Structure of 
the superconducting gap in MgB2 from point-contact spectroscopy," Supercond. Sci. Technol. 15, pp. 526- 
532, 2002; R. S. Gonnelli, D. Daghero, G. A. Ummarino, V. A. Stcpanov, J. Jun, S. M. Kazakov, and 
J. Karpinski, "Direct evidence for two-band superconductivity in MgB2 single crystals from directional 
point-contact spectroscopy in magnetic fields," Phys. Rev. Lett. 89, 247004 (2002). 

50. A. Di Chiara, F. Fontana, G. Peluso, and F. Tafuri, "Influence of the proximity effect on the conductance 
characteristics of superconducting point-contact junctions: basic assumptions," Phys. Rev. B 48, pp. 6695- 
6698, 1993; G. J. Strijkers, Y. Ji, F. Y. Yang, and C. L. Chien, "Andreev reflections at metal/superconductor 
point contacts: Measurement and analysis," ibid. 63, pp. 104510-1-6, 2001; G. Sheet, S. Mukhopadhyay, and 
P. Raychaudhuri, "Role of critical current on the point-contact Andreev reflection spectra between a normal 
metal and a superconductor," ibid. 69, pp. 134507-1-6, 2004. 

51. K. Gloos, F. Martin, C. Schank, C. Geibel, and F. Steglich, "Pressure effects in point contacts with heavy- 
fermion CeCu 2 Ge 2 ," Physica B 206 & 207, pp. 279-281, 1995. 

52. A. Golubov and F. Tafuri, "Andreev reflection in layered structures: Implications for high-T c grain-boundary 
Josephson junctions," Phys. Rev. B 62, pp. 15200-15203, 2000. 

53. Y. Kasahara, Y. Nakajima, K. Izawa, Y. Matsuda, K. Behnia, H. Shishido, R. Settai, and Y. Onuki, 
"Anomalous quasiparticle transport in the superconducting state of GeCohi^" http://arxiv.org/cond-mat/, 
pp. 0506071-1-4, 2005. 

54. N. A. Mortensen, K. Flensberg, and A. -P. Jauho, "Angle dependence of Andreev scattering at semiconduc- 
torsuperconductor interfaces," Phys. Rev. B 59, pp. 10176-10182, 1999. 

55. V. Lukic and J. Elenewski, in private communication. 

56. F. B. Anders and K. Gloos, "Towards a microscopic theory for metallic heavy-fermion point con- 
tacts," Physica B 230-232, pp. 437-440, 1997. 

57. R. Settai, H. Shishido, S. Ikcda, Y. Murakawa, M. Nakashima, D. Aoki, Y. Haga, H. Harima, and Y. 
Onuki, "Quasi-two-dimensional Fermi surfaces and the de Haasvan Alphen oscillation in both the normal 
and superconducting mixed states of CeCoIn 5 ," J. Phys.: Condens. Matter 13, pp. L627-L634, 2001. 



